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Abstract 

The production of single charged and neutral intermediate vector bosons in e + e~ 
collisions has been studied in the data collected by the DELPHI experiment at 
LEP at centre-of-mass energies between 183 and 209 GeV, corresponding to an 
integrated luminosity of about 640 pb" 1 . The measured cross-sections for the 
reactions, determined in limited kinematic regions, are in agreement with the 
Standard Model predictions. 
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1 Introduction 

The production of four-fermion final states becomes increasingly important in e + e~ 
collisions at centre-of-mass energies above the Z pole. The full set of Feynman diagrams 
must be considered, but particular topologies receive their dominant contribution from a 
subset of them (cf. Figure 1, taken from reference [1], for the standard definition of the 
different graphs). As the centre-of-mass energy increases, the dominant processes leading 
to production of vector bosons are represented by the bremsstrahlung and fusion diagrams 
with the production of a single vector boson 1 (e + e~ — > e~v e W + , e + e~ — > e + e~~f* /Z), 
and the conversion and non-Abelian annihilation diagrams leading to doubly resonant 
production (e + e~ — > W + W~, e + e~ — > ZZ). Single resonant production is dominated by 
the bremsstrahlung process, which proceeds through the scattering of a quasi-real photon 
(q^ ~ 0) radiated from an incoming e" on an e + of the other beam, i.e.: 7e + — > u e W + , 
7e + — > e + Z [2]. The resulting topology is characterized by the e~ which radiates the 
quasi-real photon being predominantly lost along the beam line. 

The integrated luminosity delivered by the LEP collider in the runs at centre-of- 
mass energies yfs = 183 — 209 GeV (LEP2) allowed, for the first time, measurements of 
the cross-section of single boson production and not just the observation of individual 
events. The evaluation of the Standard Model cross-sections for this process requires the 
computation of the full set of Feynman diagrams and, to deal with the collinear singularity 
corresponding to the electron lost along the beam line, the use of fully massive matrix 
elements. Moreover the different scales of the couplings in the process and the scale for the 
QED initial state radiation (ISR) should be properly accounted for to provide a reliable 
prediction. Therefore this process was taken as a benchmark when comparing the different 
calculations used to describe four-fermion physics at LEP2 [3]. In addition, because of 
the large missing energy in the final state, single boson production is a background when 
probing for new physics as in the search for the Higgs boson in the Hvv channel or for 
physics beyond the Standard Model [4]. Therefore the measurement of its cross-section 
is an important check that the background in these searches is correctly modelled. By 
itself, single-VT production provides access to the measurement of the trilinear gauge 
couplings at the WW'-f vertex: this measurement, in combination with other physics 
channels, has been made by the DELPHI Collaboration and is reported elsewhere [5]. 
Finally, single boson production is interesting as it will be the dominant source of weak 
boson production at a forthcoming Linear Collider. 

Single boson production is investigated in this paper in five different final states: 
e~u e qq f , e _ P e /x + i/ (U and e~v e e + v e for single-Vr production, e~e + qq and e~e + fi~fi + for 
single- Z production. Cross-sections are measured using the data collected by the DELPHI 
experiment at centre-of-mass energies ranging from 183 to 209 GeV with a correspond- 
ing integrated luminosity of about 640 pb" 1 . Compared to the results reported by the 
DELPHI Collaboration in [6], those reported here are based on a larger sample, a better 
data processing and an improved description of the simulated events (see below). There- 
fore they update and supersede those already reported in [6]. Results on single boson 
processes have been published by the other LEP experiments in [7,8] for single-W and 
in [7-9] for single- Z production. 

The criteria for the selection of the events are mainly based on the information from 
the tracking system, the calorimeters and the muon chambers of the DELPHI detector. A 
detailed description of the DELPHI apparatus and its performance can be found in [10]. 
The detector has remained essentially unchanged in the LEP2 phase, except for upgrades 

1 Charge conjugate states arc implied throughout the text. 
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of the Vertex Detector [11] and the addition of a set of scintillation counters to veto 
photons in the blind regions of the electromagnetic calorimetry at polar angles 2 9 ~ 40°, 
9 ~ 90° and 9 ~ 140°. The main tracking device was the Time Projection Chamber 
(TPC). One of the sectors (1/12 of the azimuthal acceptance) of the TPC, hereafter 
indicated as S6, was not fully operational during the last period of data taking at ^fs = 
207 GeV (about 50 pb _1 ). The data with the TPC sector down were analysed separately, 
with the performance of the analysis being evaluated on dedicated simulation samples, 
where this effect was explicitly taken into account. 

This paper is organized as follows. Section 2 introduces the definition of the single 
boson processes in terms of the four-fermion final states. In Sections 3 and 4 the se- 
lection of the events and the extraction of the single-Vr and single- Z cross-sections are 
presented. In Section 5 the comparison of measured cross-sections with the predictions 
of the Standard Model is performed for all the single boson final states. 

2 Definition of the signal and simulation 

Single boson production is investigated in this paper through four-fermion final states, 
e~P e ff and e + e~ff. These final states receive contributions not only from single reso- 
nant diagrams but also from doubly resonant production, conversion diagrams and multi- 
peripheral processes (Figure 1, according to the convention of [1]). To enhance the single 
boson production contribution, and to enable consistent comparisons and combinations 
between the experiments at LEP, the cross-sections have been defined in the limited 
kinematic regions described below. 

ev e W channel: The four-fermion final states e~u e qq' and e~v e l + vi (/ = fi, r) 

can be produced both via single-VT production, referred to as ev e W in the following, or 
via Vy-pair production. A distinctive feature of ev,W is the fact that the distribution 
of the electron direction is strongly peaked at small polar angles (9 e ) with respect to 
the incoming electron beam direction. Based on this consideration, the ev e W signal 
was defined by the complete t-channel subset 3 of the Feynman diagrams contributing 
to the e~v e qq' and e~u e l + i>i final states with additional kinematic selections to exclude 
the regions of the phase space dominated by multiperipheral diagrams, where the cross- 
section calculation is affected by large uncertainties. The signal region was therefore 
defined as follows: 

m q qi > 45 GeV/c 2 for e~u e qq', (1) 
E l+ > 20 GeV for e'v e l + Vl (l + = /i + ,r + ), 

where m, qq -i is the qq' invariant mass and Ei+ the lepton energy. Single-VT production 
accounts for more than 80% of all e~u e qq' and e~u e l + i>i events in the kinematic region 
defined above. 

ev e eh> e channel: In the kinematic region with one electron lost in the forward 

direction, this final state receives, besides single- W production, a large contribution from 
single-Z production (with Z — > u e u e ) and from the interference between single-Vr and 
single- Z processes. The contribution of this channel to the eu e W signal was also defined 

2 In the reference frame used in DELPHI the z axis was oriented along the incoming e~ beam, 9 indicated the polar 
angle and <fi the azimuthal angle. 

3 The t-channel subset consists of the bremsstrahlung, fusion and multiperipheral diagrams of Figure 1. 
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by the complete t-channel subset of the Feynman diagrams, in this case with the following 
kinematic selections: 

| cos 9 e+ 1 < 0.95, E e+ > 20 GeV and | cos e - | > 0.95. (2) 

The e~u e e + v e channel was not used in the determination of the single-Z production 
cross-section. 

Zee channel: The neutral bosons are produced in the process — > eY/Z, where 
a quasi-real photon is radiated from one of the beam electrons and scattered off the 
other beam. The signature of such events is an electron in the detector, typically of 
low energy, recoiling against the 7*/Z system, with the other electron usually lost in the 
beam-pipe. The Zee cross-section, defined for the combination of the results from all the 
LEP experiments, refers to the entire set of 48 graphs contributing at tree level to the 
e + e~ f f (/ = q, /x) final state with the following restrictions in the phase space to enhance 
the single boson contribution 4 : 

m ff > 60 GeV/c 2 and (3) 
9 e + > 168°, 60° < 6 e - < 168° and E e - > 3 GeV for a visible electron, or 
6 e - < 12°, 12° < 9 e + < 120° and E e+ > 3 GeV for a visible positron, 

the 12° (168°) being motivated by the lower angle of the acceptance for the electron 
identification of the LEP experiments. 

At yf~s = 200 GeV, within these kinematic limits, the bremsstrahlung contribution 
accounts for about 97% of e + e~qq and 67% of e + e~/i + /i~ final states. The cut on the 
invariant mass was set to 60 GeV/c 2 because it both guarantees an efficient rejection of 
the multiperipheral contribution and it provides a natural separation between the 7*ee 
and Zee regions, as it corresponds to the minimum of the differential m^j distribution. 
This paper presents, for the e + e~qq final state, besides the above defined Zee cross- 
section, a measurement of the cross-section with the same acceptance cuts for e + and e~ 
but in the invariant mass range 15 < m q q < 60 GeV/c 2 (hereafter referred to as 7*ee). 

The ev e W and Zee signal samples, as predicted by the Standard Model, were 
simulated with the WPHACT [12] event generator. For background processes, different 
generators were used: KK2f [13] for qq^), e + e~ — > /i + /i~(7) and t + t~ (7), TEEGG [14] 
and BHWIDE [15] for e+e" -> e+e"7, PYTHIA 6.143 [16] and BDK [17] for two-photon 
collisions. Fragmentation and hadronization for the KK2f and WPHACT samples were 
performed using PYTHIA 6.156. A detailed description of the simulation of four-fermion 
events at LEP2 as done in DELPHI is given in [18]. The generated signal and background 
events were passed through the detailed simulation of the DELPHI detector [10] and then 
processed with the same reconstruction and analysis programs as the data. 

3 Single- W analysis 

Both the hadronic and the leptonic final states were considered in the single-H^ anal- 
ysis. They are characterized by the presence of two hadronic jets acoplanar with the 
beam or by a single lepton with large transverse momentum, respectively. The final state 
electron is lost in the beam pipe. 

4 Diagrams involving Higgs boson exchange are neglected. Multiperipheral processes with at least one photon resolved 
are excluded from the signal definition even if they are within the accepted kinematic region. 
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3.1 Selection of hadronic events 

The experimental signature of e~v e qq' events consists of a pair of acoplanar jets. The 
undetected neutrino results in a large missing momentum at large angle to the beam 
direction. 

Other physics processes which can give rise to a similar topology are Z{j) with Z — > qq, 
WW events with at least one W decaying into hadrons, other four-fermion final states 
from neutral current processes 5 {l + l~qq, vvqq, the latter being topologically identical 
to the signal) and events induced by two-photon collisions, hereafter called two-photon 
events. Some of these processes have cross-sections larger than that of the signal by 
several orders of magnitude. Therefore the analysis was performed in two steps: a prese- 
lection of hadronic events and a final selection based on a Feed-Forward Artificial Neural 
Network [20]. 

A sample of hadronic events was preselected by requiring at least seven charged parti- 
cles to be measured in the detector. Events from Bhabha scattering were rejected by a cut 
on the total electromagnetic energy, E EM / ' yfs < 50%. The contribution from two-photon 
collisions was reduced by requiring the total visible energy to be larger than 20% of y/s 
and the total transverse energy, computed as the sum of the moduli of the momenta of 
each particle projected in the plane transverse to the beam axis, to be at least 15% of 
y/s. In addition, it was required that the half-opening angle of the cone around the beam 
axis containing 15% of the visible energy had to be larger than 10°; two-photon collision 
events are concentrated in the forward regions and have low values of this variable. The 
background from e + e~ — > qq(j) was reduced by requiring the polar angle of the missing 
momentum to satisfy the condition | cos9 miss \ < 0.98 and the acoplanarity 6 to be larger 
than 10°. Z{pl) events, with Z — > qq, were further suppressed by vetoing events with 
electromagnetic clusters with energy larger than 45 GeV or, if the ISR photon escaped 
undetected in the dead region between the barrel and end-cap electromagnetic calorime- 
ters (9 ~ 40°), by vetoing events with signals in the hermeticity counters in a cone of 
30° around the direction of the missing momentum. Selections on the maximum total 
multiplicity of charged and neutral tracks ( < 50 ) and on the visible mass (between 30 
and 100 GeV/c 2 ) were applied to suppress the residual contamination of multi-jet events 
from WW or NC processes. Finally, WW events with one W decaying to leptons were 
suppressed by requiring no identified electron or muon with energy larger than 10% or 
7.5% of y/s, respectively. Particles were identified as muons if there was at least one muon 
chamber hit associated to a track or if the size and longitudinal profile of the energy de- 
posits associated to a track in the hadronic calorimeters were consistent with a minimum 
ionizing particle. Electron identification was based on the reconstructed showers in the 
electromagnetic calorimeters associated to charged particle tracks. 

The expected composition of the residual sample after the preselection stage is shown 
in Table 1, together with the number of selected events at each centre-of-mass energy. At 
this level of the selection, the fraction of signal events is between 6% and 10% at all the 
energy points. 

The final selection of e~v e qq' events was based on a Neural Network analysis. The input 
variables were chosen to provide a good separation from the main residual backgrounds 
after preselection. The first set of variables, listed below, discriminated the signal from 
q<i{l) events: 

5 The definition of neutral current (NC) and charged current (CC) four-fermion processes of Ref. [19] are used throughout 
the text. 

6 Acoplanarity is defined as the supplement of the angle between the projections of the two hadronic jets in the plane 
transverse to the beam direction. 
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(GeV) 


e~u e qq' sij 


mal 
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NC 


qqij) 


Total MC 


Data 


183 GeV 


9.5 




66.8 


5.0 


74.2 


157.0 ±0.7 


167 


189 GeV 


32.5 




195.6 


25.7 


204.5 


462.8 ±2.0 


467 


192-202 GeV 


54.2 




276.9 


51.4 


227.3 


615.4 ± 1.7 


675 


205-207 GeV 


55.6 




240.5 


54.4 


194.9 


549.2 ±1.7 


569 



Table 1: Number of events expected from the contribution of different channels and 
observed in the data after preselection of e~u e qq' events for the different years of data 
taking. "Other CC" indicates charged current processes different from the signal, "NC" 
indicates neutral current processes. Other final states, mainly 77, contribute less than 
1% of the total number of events. The quoted errors on the total number of expected 
events ( "Total MC" ) are the ones due to limited Monte Carlo statistics. 

• effective centre-of-mass energy after ISR, Vs*, normalised to the centre-of-mass en- 
ergy, y/s [21]; 

• normalised sum of the particle momenta projected on to the thrust axis, P\ ot j y/s; 

• cosine of the polar angle of the missing momentum, | cos0 m j SS |; 

• normalised total missing momentum, P^l ss /y/s; 

• normalised total transverse momentum with respect to the beam axis, P\ ot j y/s; 

• event thrust; 

• 1 90° — 9 t hrust\, where Q thrust is the polar angle of the thrust axis. 

A second set of variables suppressed rv T qq' events, where the r lepton produces an isolated 
particle or a low multiplicity jet, and vvqq events, where the kinematic properties of the 
visible system should be consistent with the decay of a Z: 

• maximum transverse momentum of any particle with respect to the nearest jet, 
P[j ax , when the particles are clustered using the LUCLUS [22] algorithm with the 
parameter d min = 6.5 GeV/c; 

• invariant mass of the detected particles in the event rescaled by the ratio of the total 
centre-of-mass energy to the visible energy, M vis ■ y/s/E vis ; 

• Lorentz boost factor of the visible part of the event in the laboratory frame, (3 = 
P tot /E ■ 

Distributions of some of these variables, at y/s = 200 GeV, are shown in Figure 2. 
The Neural Network was trained on samples of simulated events including the signal, for 
which the output value was set to 1, and the main backgrounds, (7(7(7) and four-fermion 
processes, with output set to 0. The distribution of the Neural Network output variable is 
shown in Figure 3. The whole data sample is included in the plot. The cut on the output 
variable was set at 0.5, the value for which the product of efficiency and purity was found 
to be maximum. In the region 50 GeV/c 2 < M vis ■ y/s/E vis < 60 GeV/c 2 an excess of real 
data was found at each energy point in the final selected sample (7.4±1.3% of the total 
number of events in the real data compared to 4.2±0.2% in the simulation). To account 
for this excess, a correction of 34 fb was added to the background cross-section at each 
energy point. 

The efficiency of the selection for the signal, the expected background, the luminosity 
and the number of selected events in the data at the various centre-of-mass energies 
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Table 2: Performance (signal efficiency, e, background cross- sect ion, a^g, integrated 
luminosity, C int , and number of selected events, N data ) of the e~u e qq' event selection and 
measured cross-sections at the centre-of-mass energies considered in the analysis. The 
errors quoted on efficiencies and backgrounds are the ones due to limited Monte Carlo 
statistics. 

are reported in Table 2, together with the cross-section for the hadronic channel alone, 
evaluated from a fit of Poisson probabilities to the observed numbers of events. 

Efficiencies and backgrounds at yjs = 207 GeV were found to be compatible in the 
two periods with the sector S6 of the TPC on or off, and results have been merged in 
the table. The purity of the final selected sample is between 25% and 30%. The main 
contamination is due to WW — > ru T qq' events. 

3.2 Selection of leptonic events 

The experimental signature of the leptonic channel, e + e~ — > e~z/ e / + z//, is the presence of 
a high energy lepton accompanied by a large missing momentum and no other significant 
energy deposition in the detector. The analysis was optimised for final state leptons 
that are electrons or muons. In both channels, the contribution from ev e TV T events was 
considered as part of the background. 

The main backgrounds for the leptonic channel are the radiative production of two 
leptons, e + e~ — > l + l^( / ~f), e + e~ — > W + W~ events and two-photon collisions. 

Events were selected if exactly one well measured charged particle was reconstructed. 
The quality of the track measurement was assessed as follows: 

• relative error on the momentum, Ap/p, smaller than 100%; 

• track length greater than 20 cm; 

• polar angle 6 between 10° and 170°; 

• impact parameter in the transverse plane, \IP^\, smaller than 4 cm, and that along 
the beam direction, \IP Z \, smaller than 3 cm / sin#. 

Loose identification criteria were applied, requiring associated hits in the muon chambers 
or a significant energy deposit in the electromagnetic calorimeter. For electrons, the 
acceptance was restricted to the barrel region, | cos#| < 0.72, and the best determination 
of the electron energy was estimated by combining the momentum measurement from the 
tracking devices and the measurement of the energy deposited in the calorimeters. Any 
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Table 3: Performance (signal efficiency, e, background cross- sect ion, <Tbkg, integrated 
luminosity, £mt, and number of selected events, Ndata) of the e^v e ^i + v^ event selection at 
the centre-of-mass energies considered in the analysis. The errors quoted on efficiencies 
and backgrounds are the ones due to limited Monte Carlo statistics. 

other energy deposit in the detector unassociated to the lepton candidate was required 
not to exceed 2 GeV. In addition, the presence of tracks not fulfilling the quality criteria 
listed above was used to veto the event. The acceptance was restricted to the kinematic 
region of W decays by requiring the lepton momentum to lie below 45% of y/s and its 
transverse momentum to exceed 12% of y/s. 

A large residual contamination was still present, due to cosmic-ray events in the muon 
channel and to Compton scattering of a radiated photon in the electron channel. The 
former were suppressed by tightening the selections on the track impact parameters to 
\IPr4>\ < 0.2 cm and \IP Z \ < 2 cm for the muons. Compton scattering can mimic the 
W + — > e + u e signal when the photon balancing the electron in the transverse plane is 
lost in the dead region between the barrel and forward electromagnetic calorimeters. 
Therefore events were rejected if a signal was found in the hermeticity counters at an 
azimuthal angle larger than 90° from the electron. 

The distributions in data and simulation of the momentum of selected single muons 
and of the energy of selected single electrons are shown in Figure 4. 

The performance of the analysis at the various centre-of-mass energy values and the re- 
sults obtained are reported in Tables 3 and 4, respectively, for the e~z/ e /i + z/^ and e~v e e + v e 
channels. Due to the low statistics of selected events only the upper limits of cross-sections 
at 95% C.L. are given for each individual energy point. The limits were derived follow- 
ing a Bayesian approach from the integration of the Poissonian probabilities constructed 
with the number of events selected in the data and predicted in the simulation. For the 
electron channel a difference was found for efficiencies and backgrounds corresponding 
to the two periods at y^s = 207 GeV with TPC sector S6 on or off, and in Table 4 the 
weighted averages of the two are shown. Compatible values were found for the muon 
channel. 

3.3 Study of systematic uncertainties 

The main source of systematic uncertainty in the present measurement is the knowl- 
edge of the background level in the selected samples. In particular, as can be seen from 
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Table 4: Performance (signal efficiency, e, background cross- sect ion, <Jbkg, integrated 
luminosity, and number of selected events, Ndata) of the e~v e e + v e event selection at 
the centre-of-mass energies considered in the analysis. The errors quoted on efficiencies 
and backgrounds are the ones due to limited Monte Carlo statistics. 

Table 1, in the hadronic channel selection there is an excess of data of (5±2)% with 
respect to the expectation at the preselection stage, at which the sample consists mainly 
of background events. Rescaling the background at the final stage of the selection by 
this factor leads to an average decrease of 69±28 fb of the measured cross-section. This 
was considered as a systematic error fully correlated between the energy points. The 
correction factor added to the background cross-section, accounting for the excess of real 
data in the region 50 GeV/c 2 < M vis -y/s/ E vis < 60 GeV/c 2 , was known with a statistical 
uncertainty of ±10 fb. This leads to a systematic error on a evqq i of ±26 fb, fully correlated 
between the energy points. 

Possible inaccuracies in the modelling of background processes were evaluated by com- 
paring different Monte Carlo generators. The only effect was found in the qq{^) channel: 
using the ARIADNE [23] event generator instead of PYTHIA, the background estimate 
was 508±14 fb instead of 502±7 fb in the final e~v e qq' sample selected at 200 GeV. The 
largest of the statistical errors of the ARIADNE and PYTHIA samples was taken as a 
systematic error. The total systematic error on the background cross-section, due to the 
modelling of qq(j) process and to the limited simulation statistics (see Table 2), amounts 
approximately to ±3% in the hadronic channel. 

In the leptonic channels the systematic error on the background cross-section due to 
the limited simulation statistics (see Tables 3 and 4) amounts approximately to ±6%. 

From a comparison of dimuon events in data and simulation, the tracking efficiency, 
^tracki of DELPHI was found to be 0.5% higher in the simulation. This value was assumed 
as a systematic error. This has a negligible effect on the background, while it affects the 
selection efficiency of the signal for leptonic decays of the W. 

The uncertainty on the efficiency of the electron identification, e e , was estimated by 
comparing a sample of Bhabha events in data and simulation. The discrepancy was at 
the level of 2%. The trigger efficiency for the leptonic events is known with an error 
better than 1%. This leads to a systematic error negligible compared to the other sources 
considered. 

The luminosity is known with a total relative error of ±0.6%. 
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Table 5: Contributions to the systematic uncertainty on the e~u e qq', e~z/ e /i + z/^ and 
e~h> e e + i> e cross-sections at ^Js = 200 GeV. The selection efficiency, e, is as defined in 
tables 2, 3 and 4; e trac k is the overall tracking efficiency in DELPHI, and e e is the elec- 
tron identification efficiency. 

3.4 Total single-iy cross-section 

The total single-VF cross-section is defined as: 

Cevff = ®evqq' + 2 X O evilv + <T e uev, (4) 

where the factor two accounts for the e~v e T + v T channel, not measured in the present 
analysis, assuming /i — r universality. This assumption introduces a theoretical error at 
the level of ~ 3% on the e~v e T + v T estimation. 

The effects of the uncertainties listed in the previous section on the measurement of 
the e~V e ff cross-section at ^=200 GeV are given in Table 5. The total systematic 
error, obtained from the sum in quadrature of the individual contributions, is at the level 
of ±9%. For the measurement at the other centre-of-mass energies, the same relative 
error was assumed. 

The values of <J e vff measured at the different centre-of-mass energies together with 
their statistical and systematic errors are shown in Table 6. 

4 Single- Z analysis 

In the single 7*/Z analysis, decays of the vector boson into hadronic and final 
states were considered. Both final states are characterized by an electron scattered at 
large angle with respect to the incoming direction. The other electron, lost in the beam 
pipe, results in a missing momentum pointing along the beam direction. 

4.1 Selection of hadronic events 

The experimental signature of these events consists of a pair of jets produced in the 
hadronic decay of the ^* /Z recoiling against an electron. To maximize the sensitivity 
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Table 6: Total single-W cross- sect ion, as defined in the text (eq. 4), as measured at the 
different centre-of-mass energies considered in the analysis. The first error is statistical 
and the second one is systematic. 

of the analysis in the widest possible range of invariant masses of the ^*/Z, the event 
selection was performed in three steps: 

1. a loose preselection of events; 

2. the identification of an isolated electron; 

3. the final selection of signal events, optimized differently in two ranges of the invariant 
mass of the hadronic system, m q g, according to the most relevant background process 
in each region. 

The preselection of events consisted of the following requirements: 

• at least five charged particles in the event with at least one in the TPC with a 
measured transverse momentum larger than 2.5 GeV/c, in order to select hadronic 
events; 

• the presence of at least one electron candidate selected by requiring an energy deposit 
in the electromagnetic calorimeters E e > 3 GeV, with an associated charged particle 
with | cos6> e | < 0.985, corresponding to the acceptance of the DELPHI'S forward and 
barrel electromagnetic calorimeters; 

• in events with more than one reconstructed electromagnetic shower, the energy of the 
second most energetic shower was required to be less than 0.6E beam . This condition 
was imposed in order to reject events from Bhabha scattering. 

The electron candidates were then retained if they satisfied the following criteria: 

• in the barrel region (42° < 6 < 138°), the track parameters were required to be con- 
sistent with those of the shower measured by the electromagnetic calorimeter, with 
the additional requirement, for showers with energy higher than 30 GeV, that the 
energy deposited in the hadronic calorimeter did not exceed 10% of that deposited 
in the electromagnetic calorimeter; 

• in the forward region (10° < 9 < 32° and 148° < 6 < 170°), the energy of the shower, 
reconstructed by re-clustering the energy deposits compatible with a single electro- 
magnetic shower, was required to be compatible with the momentum of exactly one 
track measured in DELPHI'S vertex detector and very forward tracker, and with 
the momentum of no more than one track measured in the inner detector and TPC 
(see [10,11]) for detailed description of these detectors); 
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Table 7: Number of events expected from the contributions of different channels and 
observed in the data at different stages of the 7* /Zee selection (hadronic channel) for the 
different years of data taking. The column labelled "77" refers to resolved two-photon 
events. The column labelled "Others" includes other four-fermion processes, namely eeqq 
outside the signal definition and 7* /Zee with fully leptonic final state, and events from 
Bhabha scattering, the four-fermion processes supplying the more important contribution. 

• the angle, a±, of candidate tracks with respect to the closest charged particle with 
momentum p > 0.5 GeV/c had to lie in the range 15° < a± < 170°, where the upper 
limit was imposed to reject events from Bhabha scattering left in the sample at this 
stage of the selection; 

• the angle, of candidate tracks with respect to the second closest charged particle 
with momentum p > 0.5 GeV/c was required to satisfy the condition a 2 > 40°. 

Electrons from conversions or from decays were further suppressed by requiring their 
impact parameters with respect to the primary interaction vertex to be within the range 
\IPr4>\ < 0.35 cm in the transverse plane and \IP Z \ < 1 cm along the beam line. 

The charged and neutral particles were then clustered into two jets with the Durham 
algorithm [24] , excluding the tag electron. Events for which the jet resolution variable that 
separates the three-jet topology from the two-jet topology, was smaller than 10~ 4 

were rejected. For the final selection a constrained kinematic fit of the event, imposing 
energy and momentum conservation, was then performed assuming a topology of signal 
events with two jets, a visible electron and one lost along the beam line. The four- 
momentum of the invisible electron was chosen to be (0, 0, Q e E, E) with Q e the charge of 
the tagged electron. Fits with a \ 2 probability smaller than 10~ 5 were rejected. The final 
selection of signal events was then performed using the fitted kinematic variables for the 
tagged electron and the hadronic system. It was required that the following conditions 
be satisfied: 
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Table 8: Number of events expected from the contributions of different channels and 
observed in the data at the end of the 7* /Zee selection (hadronic channel) for the overall 
LEP2 sample, in the two invariant mass ranges. The column labelled "77" refers to 
resolved two-photon events. The column labelled "Others" includes other four-fermion 
processes, namely eeqq outside the signal definition and 7* /Zee with fully leptonic final 
state, and events from Bhabha scattering, the four-fermion processes supplying the more 
important contribution. 

• Q e cos9 miss > 0.95, with 9 miss being the polar angle of the missing momentum 
computed before the kinematic fit; 

• Q e cos0 e > —0.5, with 8 e being the polar angle of the tagged electron. 

The distributions of these variables after the electron identification cuts are shown in 
Figure 5 for the real and simulated data. For m qq < 60 GeV/c 2 , where the dominant 
background consisted of resolved 77 collisions, events with Q e cos9 e > 0.9 and E e > 
0.75E beam were also rejected. The numbers of selected data candidates and different 
background contributions after each selection step are shown in Table 7 for the different 
years of data-taking, while Table 8 shows the composition of the entire sample after the 
final selection in the two mass ranges. An excess of data of about 10% is observed at 
preselection level mostly due to imperfectly simulated events from Bhabha scattering. 
The efficiency of the selection on the signal, the expected background and the number of 
selected events in the data at the eight centre-of-mass energies are reported in Table 9, 
together with the evaluated cross-section. 

The distribution of the invariant mass of the hadronic system after the kinematic fit 
is shown in Figure 6 for the overall LEP2 sample. The peak in the invariant mass 
distribution around the Z mass corresponds to events for which the contribution of the 
Zee process is dominant. 

4.2 Selection of leptonic events 

The search was restricted to events with 7* jZ going into a fi + fi~ pair with invariant 
mass above 60 GeV/c 2 . The general features are the same as for the hadronic channel 
with jets replaced by muons. Thus a three-track signature, of two high-momentum muons 
and one e + or e~ , scattered at large angle, is expected in the detector. The signal selection 
criteria on angular distributions were similar to those used in the hadronic channel. 

In the preselection the event was required to have exactly three tracks fulfilling the 
following criteria: 

• fractional error on the momentum Ap/p < 50%; 

• impact parameter in the transverse plane \IPr^\ < 0.5 cm and along the beam 
direction \IP Z \ < 3 cm; 

• at least one associated hit in the Vertex Detector. 
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0.161&M ± 0.02 


202 


29.0±0.8 


0.026±0.005 


40.3 


2 


0.08+^ ± 0.03 


205 


29.8±0.7 


0.021±0.004 


75.9 


12 


± 0.03 


207 TPC OK 
TPC S6-off 


28.3±0.9 
27.4±0.7 


0.019±0.004 
0.030±0.014 


84.1 
51.4 


4 
9 


0.25l£o9 ± 0.02 


Y/Z -> qq (m q q > 60 GeV/c 2 ) 


(GeV) 


(%) 


O'bkg 
( P b) 


(pb" 1 ) 


Ndata 


(pb) 


183 


27.2±0.4 


0.078±0.017 


52.0 


12 


0.56t^ ± 0.07 


189 


27.8±0.4 


0.068±0.007 


153.5 


38 


0.64+^ ± 0.04 


192 


28.1±0.4 


0.063±0.006 


25.1 


6 


0.63t&$ ± °- 04 


196 


28.8±0.3 


0.060±0.006 


75.9 


19 


0.66t^ ± 0.04 


200 


29.7±0.5 


0.072±0.006 


82.8 


20 


0.571^ ± 0.04 


202 


30.5±0.4 


0.066±0.006 


40.3 


5 


0.191^ ±0.02 


205 


30.7±0.3 


0.072±0.006 


75.9 


14 


0.371^ ±0.03 


207 TPC OK 
TPC S6-off 


31.0±0.3 
29.4±0.3 


0.068±0.006 
0.060±0.004 


84.1 
51.4 


22 
15 


omt° ii ± 0.03 


Y/Z -> (m M+/1 - > 60 GeV/c 2 ) 


(GeV) 


(%) 


®bkg 
(fb) 


r 

*- / int 
(fb" 1 ) 


Ndata 


a 

(fb) 


183 


27.4±1.1 


0.6±0.2 


0.0540 


1 


< 319 at 95% C.L. 


189 


26.2±1.0 


1.1±0.4 


0.1581 


5 


C 250 at 95% C.L. 


192 


26.3±1.0 


0.7±0.2 


0.0258 





; 441 at 95% C.L. 


196 


26.7±1.0 


1.3±0.3 


0.0769 


2 


< 301 at 95% C.L. 


200 


27.2±1.0 


1.2±0.3 


0.0843 


1 


< 203 at 95% C.L. 


202 


26.7±1.0 


0.9±0.2 


0.0411 





< 273 at 95% C.L. 


205 


26.4±1.0 


0.6±0.2 


0.0767 


1 


< 232 at 95% C.L. 


207 TPC OK 
TPC S6-off 


26.1±1.0 
27.5±1.0 


1.1±0.3 
0.8±0.2 


0.0874 
0.0544 


1 
2 


< 201 at 95% C.L. 



Table 9: Performance (signal efficiency, e, background cross-section, a bkg , integrated 
luminosity, Ci nt , and number of selected events, N^ata) of the 7* /Zee event selection at 
the centre-of-mass energies considered in the analysis. The period with TPC sector 6 
down is indicated as "207 TPC-S6 off". Cross-sections for Y/Z — > are expressed 

in femtobarns due to smaller values. 
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7*/Zee (e + e )bk g Others 


Total MC 


Data 


183 GeV 








Preselection 


0.87±0.03 6.3±0.1 0.6±0.1 


7.8±0.1 


4 


Final Selection 


0.60±0.02 0.008±0.003 0.02±0.01 


0.63±0.03 


1 


189 GeV 








Preselection 


2.6±0.1 19.2±0.4 1.8±0.2 


23.6±0.4 


24 


Final Selection 


1.7±0.1 0.03±0.01 0.14±0.06 


1.91±0.09 


5 


192-202 GeV 








Preselection 


3.8±0.1 26.3±0.3 2.4±0.1 


32.5±0.3 


21 


Final Selection 


2.7±0.1 0.04±0.01 0.21±0.03 


2.93±0.06 


3 


2000 205-207 GeV 








Preselection 


3.7±0.1 23.8±0.3 2.3±0.1 


29.8±0.3 


39 


Final Selection 


2.6±0.1 0.04±0.01 0.15±0.03 


2.79±0.06 


4 



Table 10: Number of events expected from the contributions of different channels and 
observed in the data at different stages of the 7*/Zee selection (leptonic channel) for the 
different years of data taking. The column labelled "(e + e _ fi + ^~)bk g n shows the numbers 
for the background events coming from all processes with e + e~fi + fi~ in the final state not 
fulfilling the signal definition criteria. All other background sources are collected inside 
the column labelled "Others". 

The sum of the charges of the three particles was required to be ±1. Possible photon 
conversions were removed according to the procedure described in [10] and requiring in 
addition the minimum opening angle of any track pair to be larger than 5°. 

Since the event topology is simple, the particle identification required at least two 
tracks to be identified as leptons (/i or e) and at least one of them to be a muon. For 
muon identification loose criteria were applied, as in the case of single-VF production 
(see Section 3.2). The flavour of the possible unidentified track was inferred from partial 
information taking into account the combination of the charges of the observed particles. 
In the case of n + x~e ± or x + /i~e ± , the unidentified track x was treated as \i. For \i + n~x ± 
the track x was taken as e ± . Since the efficiency of the identification of the electrons was 
smaller than for the muons, a majority of events with the n + fi~ pair detected and an 
unidentified electron was accepted this way. This reduced dramatically the sensitivity of 
event selection, and hence the loss of efficiency, to the efficiency of electron identification 
(less than 5% drop of signal selection efficiency was observed after forcing the electron 
track to be always unidentified). Events with two tracks identified as electrons were 
rejected. At the preselection stage, the momentum of the e ± candidate had to be greater 
than 2 GeV/c, and the invariant mass of the fi + fi~ pair greater than 20 GeV/c 2 . 

Due to the stringent cut on low multiplicity of the event, the data reduction factor 
was large. For all energy points, 88 events were preselected and 94.0±0.6 events were 
expected. At this stage most of the events came from other neutral current four-fermion 
processes with e + e~fi + fi~ in the final state but outside the kinematic limits of the signal 
definition (see Table 10). The remaining contributions came mainly from the neutral 
current four-fermion processes (e + e~ — > l^l^l^l^ excluding the e + e~[i + fi~ case), from 
two-fermion processes (e + e~ — > and e + e~ — > t + t~{j)) and a small fraction from 

e+e" -> W + W-. 
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A kinematic fit was performed before applying the final selection cuts to the data. An 
electron lost along the beam line and no missing momentum in the transverse plane were 
assumed. The invariant mass of the pair was recalculated if the probability of the 

kinematic fit was above 0.001. The high purity of this channel implies that events with 
lower fit probability still correspond to e + e~fi + fi~ production, but with a poorly measured 
visible electron; such events were therefore retained and the original uncorrected 
invariant mass was kept. In agreement with the signal definition the ji + ji~ invariant mass 
was then required to be greater than 60 GeV/c 2 , the momentum of the observed electron 
greater than 3 GeV/c and its polar angle 9 e to satisfy the condition Q e cos9 e > —0.5, 
with Q e being the charge of the observed electron. 

Finally the allowed angular ranges for the direction of the Z/Y momentum and missing 
momentum were defined by the following conditions: 

• Q e cos9 fM +^- < —0.8, with 9^+^- being the polar angle of the fi + fi~ system; 

• Q e cos9 miss > 0.8, with 9 miss being the polar angle of the missing momentum com- 
puted before the kinematic fit. 

After the final selection the background contribution is expected to be less than 10% of 
the total selected events. The remaining background from processes with e + e~/i + /i~ in 
the final state which was dominant at the preselection level was reduced to about 1%. 

The efficiency of the selection of the signal, the expected background and the number 
of selected events in the data for all centre-of-mass energies are reported in Table 9. Due 
to the low statistics of selected events only the upper limits of cross-sections at 95% 
C.L. are given for each individual energy point. The limits were derived following a 
Bayesian approach from the integration of the Poissonian probabilities constructed with 
the number of events selected in the data and predicted in the simulation. In total 13 
events were selected and 8.3±0.1 events were expected from data in the energy range 
from 183 GeV to 207 GeV. The n + fi~ invariant mass distribution is shown in Figure 6. 

The distributions of the energy and of the signed angle, Q e cos9 e , of the tag electron 
after the kinematic fit for hadronic and ji~ events with nifj > 60 GeV/c 2 , are shown 
in Figure 7 summed over all the centre-of-mass energies. The observed spectra are in 
agreement with the predictions from the simulation. 

4.3 Systematic uncertainties 

The measurement uncertainty is dominated by the limited real data statistics. 

In the hadronic channel three sources of systematic errors were considered: the effi- 
ciency in the electron selection procedure, the modelling of the contribution from two- 
photon events, which represents the largest background component in the low invariant 
mass region, and the modelling of the fragmentation in the qq("i) contribution, which 
represents the largest background component in the high invariant mass region. 

The uncertainty on the efficiency of the electron identification was estimated by com- 
paring the number of selected events in the data and in the simulation for a sample 
enriched in WW events (about 85% purity) with at least one of the two W's decaying, 
directly or in cascade, into a final state containing an electron. The same criteria for 
electron identification and isolation were adopted as in the Zee analysis, but the signal 
selection criteria were changed to maximize the product efficiency times purity of the WW 
selection. Assuming the Standard Model prediction for the WW cross-section, includ- 
ing the O(a) electroweak corrections via the so-called Leading Pole Approximation [3], 
and attributing the entire discrepancy between the observed and the expected number of 
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events to the different electron identification efficiency in the data and in the simulation, 
the relative difference in the efficiency was found to be Ae e /e e = (—2.2 ±3.6)% where the 
error accounts both for the data and the simulation statistics. The error on the difference 
was used for the computation of the systematic error. 

The uncertainty in modelling of two-photon events could arise from the bad modelling 
either of the direct or of the resolved photon contribution. As described in [18], in the 
region of single tag the direct component was simulated using the WPHACT generator 
and the resolved component using PYTHIA 6.143. To match the direct and resolved 
components in the region m qq < 40 GeV/c 2 the WPHACT generator was run with 
constituent quark masses. The direct component of single tag events with m qq > 40 
GeV/c 2 was instead simulated with the WPHACT generator using current algebra quark 
masses. To gauge the effect of the different quark masses for the single tag low mass 
direct component, a fully simulated sample with current algebra quark masses only was 
used to evaluate the effect both on signal efficiency and on the background cross-section. 
The change in the quark mass does not affect the 7* /Zee signal at any stage of the 
selection, while the background at the end of the selection is increased by about 5 fb in 
each invariant mass region. Concerning the resolved photon component, at y/s = 200 
GeV the cross-section of this background amounts, at the generator level, to about 10 fb 
in the 7*ee signal region and 17 fb in the Zee one. A different generator, TWOGAM [25], 
predicts background cross-sections of 5 fb and 9 fb, respectively, in these two regions. As 
the topologies of the resolved photon and 7* /Zee events are similar, the same selection 
efficiency was assumed for this background and the signal. Therefore the difference 
between the PYTHIA and TWOGAM predictions at the generator level, which is stable 
in the range y/s = 183 — 207 GeV, was taken as systematic error on the measured cross- 
section. 

The uncertainty in modelling the fragmentation and hadronization in qq{^f) events 
was evaluated using a simulation sample produced with the ARIADNE generator. The 
background cross-sections were found to be larger, but within the statistical error, leading 
to a decrease of the measured cross-sections of 3±5 fb in the low invariant mass region 
and 7±13 fb in the high invariant mass one. The largest of the statistical errors from the 
ARIADNE and PYTHIA samples, 3 fb in the 7*ee signal region and 9 fb in the Zee one, 
were taken as systematic error. 

These three systematic uncertainties, together with the error on the luminosity, were 
taken as fully correlated at the different centre-of-mass energies, while the errors on 
the background cross-section and on the signal efficiency due to the limited simulation 
statistics were considered uncorrelated among the different energies. 

The contributions of the sources of systematic uncertainty in the hadronic channel 
at 189 GeV are summarized in Table 11. The total systematic uncertainty amounts to 
±10% in the region 15 < m qq < 60 GeV/c 2 and to ±6% for m qq > 60 GeV/c 2 . 

The contributions of the different sources of systematic errors in the leptonic channel 
are summarized in Table 12. The main source of systematic error is the limited simulation 
statistics, both for the signal and for the background. The effect of the uncertainty on the 
efficiency of the electron identification was measured to be negligible using relaxed iden- 
tification criteria. The total systematic uncertainty amounts to about ±5% per energy 
point. Assuming no energy correlation of the systematic errors, the overall systematic 
uncertainty on the energy averaged cross-section was estimated to be ±2.5%, an order of 
magnitude smaller than the statistical uncertainty. 
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Svstpmatir pffprt 


Error on a (pb) 




1 ^ <r m - <r (SD HpV /V 2 

ItJ \ ''Lqq \ UU vjrc V /C 


m - fin ripV /V 2 


Ae e 


0.009 


0.022 


Aa bkg (77) direct 


0.005 


0.005 


Aa bk g (77) resolved 


0.005 


0.008 


A/Jbkg (997) fragmentation 


0.003 


0.009 


Ae due to simulation statistics 


0.007 


0.008 


AfTbfcc, due to simulation statistics 


0.020 


0.026 


Luminosity ±0.6% 


0.002 


0.005 


Total 


0.024 


0.038 



Table 11: Contributions to the systematic uncertainty on the 7* /Zee cross-sections in 
the hadronic channel, in the two ranges of invariant mass of the hadronic system, at 
y/s = 189 GeV. The selection efficiency, e, is as defined in table 9 and e e is the electron 
identification efficiency. 



Systematic effect 


Error on (fb) 


Ae due to simulation statistics 
Aa bkg due to simulation statistics 


1.8 
1.4 


Luminosity ±0.6% 


0.7 


Total 


2.4 



Table 12: Contributions to the systematic uncertainty at ^fs = 189 GeV on the predicted 
7*/Zee cross-section in the leptonic channel. The selection efficiency, e, is as defined in 
table 9. The systematic errors were conservatively considered to be the same for all cen- 
tre-of-mass energies. The systematic due to the uncertainty on the electron identification 
efficiency was measured to be negligible. 

5 Combined single boson cross-sections 

The measured values for single boson cross-sections are compared with the Standard 
Model predictions obtained with WPHACT [12] as a function of the centre-of-mass en- 
ergy. This dependency is shown in Figure 8 for single-VT and single-^ production. The 
theoretical uncertainty on the predictions amounts to 5%. The overall compatibility 
with the Standard Model was checked by considering the ratio R of the measured to 
the predicted cross-sections. At each energy point a Poissonian probability function was 
constructed based on the number of observed events, the number of expected background 
events and the signal extraction efficiency. A maximum likelihood fit to the data of the 
global probability function, being the product over all probability functions for individ- 
ual energies convoluted with a multidimensional Gaussian describing the correlated and 
uncorrelated systematic errors, was performed. The results were: 
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R(eis e qq') = 1.20 ± 0.18 (stat.) ± 0.14 (syst.), 
R[ev e iiv^) = 1.06±8:1I (stat.) ±0.03 (syst.), 
R[ev e ev e ) = 1.07tg;l (stat.) ± 0.09 (syst.), 
R(eeqq) = 1.22lg;l| (stat.) ± 0.06 (syst.) 15 < m qq < 60 GeV/c 2 , 

R(eeqq) = l.OOtan ( stat ± °- 05 ( s y st m m > 60 GeV/c 2 , 

R(eefxfL) = 1.59t° H ( sta t-) ± 0.03 (syst.) m M+M - > 60 GeV/c 2 , 

where the systematic error represents only the experimental contribution. The values 
found show a good agreement with the Standard Model predictions. 
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Abelian Classes 




Fusion 



Figure 1: Four-fermion production classes of 
convention of [1]: B = Z,j and Bi,B 2 , B 3 = 
exchange are not shown. 




Annihilation 



diagrams in e + e annihilation following the 
Z, 7, W ± . Diagrams involving Higgs boson 
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Figure 2: ev e W channel (W — > qq') at ^/s = 200 GeV. Distribution of some Neural Net- 
work input variables, as defined in the text, in the real data (points with error bars) and 
in the simulation of the Standard Model predictions (histograms) after the preselection 
stage (see text). The distributions of these variables for the ev e W signal are shown as well 
(filled histograms). The | cos# m i SS | distribution (top-right) stops at 0.98, as explained in 
the text. 
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DELPHI 




NN output 



Figure 3: ev e W channel (W — > g<f) summed over all centre-of-mass energies: distribution 
of the Neural Network output variable in the real data (points with error bars) and in 
the simulation of the Standard Model predictions (histograms). The filled histogram 
represents the single-V^ signal, the open area is the background expectation. The arrow 
indicates the cut applied on this variable for the final event selection. 
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Figure 4: ev e W channel (W — > summed over all centre-of-mass energies: momentum 
distribution of the muon (top) and energy distribution of the electron (bottom) in the real 
data (points with error bars) and in the simulation of the Standard Model predictions 
(histograms) for the events selected at the end of the analysis. The filled histograms 
represent the single-W signal, the open area is the background expectation. 
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Figure 5: 7*/Zee channel (j*/Z — * qq) summed over all centre-of-mass energies: distri- 
butions of the variables used for the signal definition at the reconstruction level after the 
"electron identification" step (see Section 4.1), in the real data (points with error bars) 
and in the simulation of the Standard Model predictions (histograms). The 7*/Zee signal 
is defined in the kinematic region described in Section 2. 
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Figure 6: 7* /Zee channel summed over all centre-of-mass energies: invariant mass distri- 
bution of the ^*/Z system in the real data (points with error bars) and in the simulation 
of the Standard Model predictions (histograms) for hadronic (top) and fi + fi~ (bottom) 
final states, in the selected signal sample. 
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Figure 7: 7*/Zee channel summed over all centre-of-mass energies: energy spectrum 
(top) and signed angle Q e cos 9 e , (bottom) of the tag electron for hadronic and final 
states with m^j > 60 GeV/c 2 , in the selected signal sample. The points with error bars 
represent real data, the histograms the simulation. 
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Figure 8: Cross-sections as a function of y/s for the ev e W channel (top) and the Zee 
channel Zee — > e + e~qq (bottom). The solid curves are the Standard Model predictions, 
with a 5% uncertainty band, computed with WPHACT [12]. 



